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1.0  RESEARCH  OBJECTIVES 


We  had  two  objectives  for  this  effort.  First,  we  wanted  to  develop  a  systematic  approach 
to  determine  radiation  sensitivity  in  organic  semiconductor  based  devices,  using  photo-voltaic 
devices  as  a  matter  of  expediency.  Our  second  objective  was  to  apply  the  test  and  modeling 
techniques  to  characterize  an  array  of  organic  based  devices  (FETs,  photo-cells,  etc.)  developed 
through  other  projects. 

2.0  SUMMARY  OF  FY11  PROGRESS 

In  FY 1 1  we  made  significant  advances  in  the  understanding  and  modeling  of  radiation 
effects  in  organic  semiconductor  devices.  We  used  device  modeling  to  simulate  the  behavior  of 
P3HT:PCBM  based  photo-cells,  we  developed  a  qualitative  description  of  the  open-circuit 
voltage  (Voc)  variations  and  photo-induced  carrier  relaxation/recombination  time  constant  (x) 
variations  with  incident  light  level,  and  we  approximated  the  effects  of  irradiation  by  a  build-up 
of  charge  near  the  anode  contact  leads.  This  last  approximation  provides  a  satisfactory 
description  of  the  experimentally  observed  behavior  of  Voc  and  x. 

3.0  BACKGROUND 

It  remains  our  conviction  that  organic-semiconductor-based  photo-cells  may  be  usefully 
employed  as  power  sources  in  space,  particularly  for  short-term,  rapidly-assembled,  low  earth 
orbit  (LEO)  missions.  Though  they  have  disadvantages  over  classic  inorganic-based  cells  in  that 
their  photo-conversion  efficiency  is  low  (less  than  8%  presently),  their  specific  power  (power  per 
unit  mass)  can  be  substantially  higher  given  the  thin  film,  low  mass  nature  of  the  cells.  One  can 
easily  imagine  their  incorporation  on  the  body  of  a  satellite  or  on  the  surface  of  an  inflatable, 
easily  deployed  structure.  However,  there  have  been  few  studies  of  the  effects  of  radiation 
(primarily  the  ionizing  component  due  to  high  energy  protons  and  electrons)  to  which  photo¬ 
cells  are  subjected  in  a  LEO.  For  this  effort  we  used  an  x-ray  source  to  simulate  the  ionizing 
radiation  in  space  both  to  characterize  the  induced  damage  and  to  serve  as  a  basis  for  theoretical 
calculations  of  the  physics  of  radiation  damage.  Therefore,  this  report  is  divided  into 
experimental  and  theoretical  sections,  albeit  overlapping. 

4.0  EXPERIMENT 

Early  in  this  effort  we  introduced  the  short  light  pulse  method1  in  which  a  light  pulse  of 
approximately  2  ps  duration  is  superimposed  upon  a  constant  but  adjustable  background  light 
level.  That  background  light  is  used  to  set  the  open-circuit  voltage  (Voc)  level.  The  short  light 
pulse  causes  a  small  perturbation  in  photo-induced  carrier  density  (An),  which  then  induces  a 
variation  in  Voc  (AVoc).  That  variation  subsequently  relaxes  in  a  time,  x,  characteristic  of  the 
recombination  of  the  carriers.  We  measured  Voc  using  a  high  input  impedance  voltmeter  and 
AV0c(t)  using  a  high  input  impedance,  digitizing  oscilloscope.  This  method  to  determine  x  and 
Voc  as  a  function  of  light  intensity  or  radiation  dose  was  employed  throughout  the  experiments 
reported  here. 
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Electron-hole  pairs  are  generated  when  the  cell  is  illuminated  by  a  light  source  whose 
photon  energy  is  greater  than  the  band-gap  energy.  The  electron-hole  pair  density  increases  with 
illumination  time  until  the  recombination  rate  equals  the  generation  rate.  In  our  experiments,  the 
excess  carriers  were  generated  by  a  short  light  pulse  resulting  in  an  exponential  generation  and 
decay  of  excess  carriers.  While  exposing  the  cells  to  these  light  pulses  we  measured  the  open 
circuit  voltage  as  a  function  of  time  and  used  the  Voc  rise  and  fall  times  as  a  measure  of 
t.  During  these  experiments  we  observed  that  x  changes  with  excess  carrier  density — an 
observation  that  will  be  important  in  the  following  discussions. 

Figure  1  shows  the  results  of  four  experiments.  For  the  first  experiment  we  varied  Voc  by 
adjusting  the  background  light  level  while  measuring  x.  In  the  other  three  experiments  Voc  was 
set  to  a  specific  value  (0.34,  0.40  and  0.44  V)  and  x  was  measured  while  the  radiation  was 
increased  from  zero  to  300  krad  (SiCfi).  For  these  last  three  experiments  x  appeared  to  be 
independent  of  dose.  This  result  led  us  to  conclude,  and  subsequently  confirm,  that  there  is  a 
build-up  of  charge  with  radiation  dose  near  the  PEDOT  PSS  (Poly(3,4-ethylenedioxythiophene) 
poly(styrenesulfonate))  anode  contact.  And,  because  x  does  not  change  with  radiation  dose  we 
assume  that  carrier  recombination  is  not  influenced  by  this  charge  build-up. 


Figure  1.  The  effect  of  radiation  on  photo-induced  carrier  relaxation/recombination  time 

constant  (t). 


4.1  New  Radiation  Data 

In  view  of  the  data  shown  in  Fig.  1  we  proceeded  to  adjust  different  variables  such  as  the 
P3HT:PCBM  ratio,  then  studied  the  effect  of  radiation  for  accumulated  doses  up  to  50  krad 
(SiCk)  as  shown  in  Fig.  2.  The  initial  value  of  Voc  at  dose  =  0  krad  (SiC>2)  was  the  maximum 
value  attainable  for  that  sample  with  the  available  halogen  light  source.  Strikingly,  high  weight- 
ratio  samples  appeared  relatively  radiation  insensitive  while  those  optimized  for  photo¬ 
conversion  efficiency  unfortunately  showed  a  maximum  for  the  ratios  measured. 


2 

Approved  for  public  release;  distribution  unlimited. 


0.46 
0.45 
0.44 
0.43 

>  0.42 
o 
o 

>  0.41 
0.40 
0.39 
0.38 

0  10  20  30  40  50 

Accumulated  Dose  (krad  (SiC>2)) 


■  G11  (1:1) 

•  G1  (10:1) 

A  G9  (2:1) 

▼ 

▼  G4  (4:1) 

▲ 

G15  (1:2) 

•  • 


Figure  2.  Effect  of  accumulated  radiation  on  Voc  for  different  weight  ratios  of  P3HT  to 

PCBM. 


Note  that  in  Fig.  3  that  the  AV0C  values  are  shown  for  a  total  accumulated  dose  of  50  krad 
(SiCh)  and  that,  since  the  data  in  Fig.  2  corresponds  to  essentially  straight  line  plots,  one  can 
determine  the  slope,  AVoc/ADose  by  dividing  the  experimental  data  in  Fig.  3  by  50  krad  (SiCF). 


Weight  Ratio  (P3HT:PCBM) 

Figure  3.  AV0C  plotted  as  a  function  of  weight  ratio. 


We  previously  concerned  ourselves  with  the  potential  issue  of  radiation  induced  charging 
of  the  glass  cover  slide  and  its  possible  role  in  inducing  “effective”  shifts  in  the  Voc  parameter  as 
shown  in  Fig.  1.  Two  facts  now  emerge  to  counter  this  possibility.  First,  the  data  in  Fig.  2  is 
inconsistent  with  such  a  charging  phenomena  since  all  of  the  curves,  irrespective  of 
P3HT:PCBM  weight  ratio,  would  be  expected  to  show  some  dose-dependent  variation  of  Voc  and 
this  is  not  the  case.  Furthermore  we  have  studied  dose  dependence  of  Voc  in  1:1  by  weight 
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samples  having  cover  slides  differing  in  thickness  such  that  their  x-ray  attenuation  factors 
differed  by  a  factor  approximately  4.  For  the  same  accumulated  x-ray  dose,  as  received  by  the 
photo-cell  post  attenuation,  we  observed  the  same  Voc  variation.  This  eliminates  the  possibility 
of  charging  effects  coming  from  x-ray  absorption  in  the  glass. 


4.2  Dose  Rate  Effect 

One  of  the  issues  that  was  not  clear  to  us  was  whether  or  not  different  dose  rates  had  any 
effect  on  cell  charging  since,  as  shown  in  Fig.  4,  we  have  evidence  for  a  relaxation  phenomenon 
when  irradiation  is  halted  (times  after  ~  50  minutes  in  Fig.  4).  We  will  demonstrate  that  the 
presence  of  a  relaxation  term  necessarily  results  in  dose-rate  dependence.  Currently,  Voc  is  the 
observable  parameter  in  Fig.  4,  but  we  are  more  interested  in  the  charging  of  the  device  resulting 
from  irradiation.  As  an  initial  starting  point,  we  assumed  a  first  order  charge  generation 
mechanism  (the  first  term  on  the  right  side  of  Eq.  1)  and  a  simple  linear  relaxation  term.  This 
relation  is  consistent  with  indirect  recombination  in  the  limit  of  low-level  injection  . 
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Figure  4.  Voc  versus  irradiation  time  up  to  a  cumulated  incident  dose  ~  500  krad(Si02). 
Subsequently  the  device  was  allowed  to  relax  (t  >  50  min). 

TO 

This  yields  the  following  solution  for  the  evolution  of  the  charge  density  as  a  function  of  time: 

<2> 


Relaxation 
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In  this  last  equation  P(t)  is  the  trapped  charge  concentration,  D  is  the  dose  rate,  ap  is  the  charge 
capture  cross  section  (neutral  traps),  cj)  is  the  electron-hole  escape  probability,  P0  is  the  maximum 
density  of  traps,  and  k  is  a  relaxation  parameter.  Equation  2  suggests  that  there  is  dose-rate 
dependence  through  two  factors:  a)  if  k  >  Dapc()  or  k  ~  Dapc(),  then  the  saturation  level  is  less 
than  P0  and  is  reduced  to  Dap(j)P0/k  for  those  instances  in  which  k  »  Dap(j)  and  b)  the  initial 
slope  (dP/dt  |t — >0)  is  equal  to  Dap(t)Po.  To  investigate  this,  different  dose  rates  were  used  with 
1:1  devices.  Figure  5  shows  the  relationship  between  Voc  and  accumulated  dose  for  different 
dose  rates  in  these  1 : 1  samples.  This  analysis  assumes  linearity  and  constancy  of  coefficients  in 
the  linear  ordinary  homogeneous  differential  equation.  We  approximated  AVoc  as  AP,  although 
this  assumes  a  constant  coefficient  linear  relationship  between  charge  density  and  Voc  which  is 
not  necessarily  correct.  This  point  will  be  addressed  later. 
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Figure  5.  AV0C  plotted  as  a  function  of  accumulated  dose  for  1:1  samples.  Attenuation 
factor  is  being  accounted  for  in  the  dose  rates.  All  dose  rates  are  in  krad/min  (Si02). 

From  Fig.  5  it  would  appear  that  there  is  some  dose-rate  effect  on  AV0C — although  not  a 
linear  function  of  dose  rate  as  we  might  have  anticipated.  However,  it  must  be  noted  that  we 
were  very  limited  in  terms  of  the  lowest  dose  rate  we  could  use.  For  example,  the  data  taken 
using  a  dose  rate  of  0.49  krad  (SiC>2)  required  the  x-ray  machine  to  function  for  more  than  100 
minutes,  which  is  not  desirable.  Using  a  higher  dose  rate  did  not  produce  abnormal  effects  on 
Voc  and  the  results  were  about  the  same  as  when  the  dose  rate  was  7.3  krad/min  (SiCF). 
However,  we  believe  that  the  change  in  Voc  will  reach  a  saturation  limit  at  higher  dose  rates. 
This  could  explain  why  the  change  in  Voc  was  slightly  higher  for  the  3.3  krad/min  (SiCF)  dose 
rate  than  for  the  6.5  krad/min  (SiC>2)  dose  rate.  Another  cause  could  be  noise  in  the  experiment. 

4.3  Relationship  between  Charge  and  Accumulated  Dose 

As  mentioned  previously,  we  plotted  AV0C  as  a  function  of  dose  in  Fig.  5  but  would  have 
preferred  to  plot  charge  density  as  a  function  of  dose.  Using  theoretical  modeling  (to  be 
discussed)  we  can  translate  Voc  into  a  charge  density  and  re-plot  Fig.  5  as  Fig.  6.  A  linear 
relationship  between  charge  density  and  accumulated  dose  appears,  showing  that  the 
approximate  form  used  to  generate  Fig.  5  was  adequate. 
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Figure  6.  Charge  plotted  as  a  function  of  accumulated  dose. 


4.4  Ancillary  interface  modification  studies 

We  previously  attributed  variations  in  Voc  in  1:1  samples  as  being  due  to  charge  buildup 
at  the  anode  interface.  We  also  carried  out  experiments  on  1:1  samples  in  which  the  cathode 
interface  was  modified  by  addition  of  thin  films  of  ZnO,  LiF,  TiCfi/Cs,  and  standard  Ca/Al  or 
pure  Al.  These  cells  were  again  subjected  to  varying  levels  of  optical  irradiation  from  a  halogen 
lamp  and  the  voltage  relaxation  time  was  measured  versus  Voc  (see  Fig.  7). 
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Figure  7.  Experimental  results  of  the  relaxation  time  measurements  (logarithmic  scale)  of 
1:1  P3HT:  PCBM  samples  with  different  interfaces  at  the  cathode. 


Following  Koster,  et  al.,  one  can  express  x  in  the  fonn: 


gVpc 
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(3) 


where  q  is  the  elementary  charge,  k  the  Boltzmann  constant,  T  the  temperature,  Egap  is  the  energy 
difference  between  the  highest  occupied  molecular  orbital  (HOMO)  of  the  electron  donor  and  the 
lowest  unoccupied  molecular  (LUMO)  of  the  electron  acceptor,  y  the  recombination  constant, 
and  Nc  the  density  of  states  in  the  conduction  band.  From  this  expression  we  anticipate  a  linear 
relationship  between  ln(x)  and  Voc  with  a  slope  of  -  q/kT  and  this  should  be  independent  of  the 
interface.  This  is  clearly  not  the  case,  indicating  that  modification  of  the  cathode  interface  varies 
the  exponential  factor  which  should,  in  fact,  be  a  constant.  At  the  present  time  we  have  no 
explanation  for  this  phenomenon  but  hope  that  theoretical  modeling  will  provide  some  insight. 
Irradiation  studies  of  these  interface-modified  systems  should  also  help  explain  the  phenomenon. 

5.0  MODELLING 

In  the  following  sections  we  describe  our  efforts  to  model  the  effects  observed 
experimentally. 

5.1  Simulations  (General  Information) 

To  simulate  the  organic  photocells  we  used  device  simulation  software  designed  for 
conventional  semiconductors  such  as  silicon  and  gallium-arsenide  (GaAs).  Two  codes  were 
used:  Radiation  Effects  in  Oxides  and  Semiconductors  (REOS)4  (developed  at  Sandia  National 
Laboratories)  and  Silvaco  Atlas5  which  is  a  commercial  device  simulation  code.  Silvaco  Atlas 
works  by  solving  the  drift-diffusion  equation  and  Poisson  equations  simultaneously  in  order  to 
predict  the  carrier  densities  and  potentials  within  the  device  under  known  biasing.  It  then  uses 
those  values  to  calculate  the  current. 

For  our  initial  approach  we  used  REOS  to  model  the  organic  solar  cell  as  a  layered  PN 
junction  similar  to  a  silicon  solar  cell  but  with  changed  physical  parameters  to  represent  the 
lower  mobilities  and  different  band  gap  of  the  materials.  This  approach,  however,  is  physically 
incorrect  since  the  blended  organic  material  is  not  two  distinct  layers  with  high  doping.  Later  we 
switched  to  Silvaco  based  on  an  already  implemented  Langevin  recombination  model  and  the 
advantage  of  highly-debugged  commercial  code. 

Presently,  the  device  is  modeled  as  a  photovoltaic  device  that  is  200  nrn  thick,  with  an 
effective  area  of  2  mm  by  1  mm  or  2  mm  .  To  accomplish  the  simulation  in  Silvaco,  the  device 
structure  is  modeled  as  a  continuous  device  with  uniform  cross  section  from  the  cathode,  through 
the  PCBM:P3HT  (also  termed  the  blend)  to  the  PEDOT,  (which  functions  as  the  anode)  with  no 
variations  in  the  blend  or  the  cross  sectional  area.  The  blend  is  treated  as  a  single  material  with  a 
single  set  of  physical  parameters  in  which  a  built-in  potential  is  created  by  the  work  functions  of 
the  two  contacts — similar  to  the  method  used  in  Reference  6.  The  cathode  is  presently  assumed 
to  be  aluminum  and  the  effective  work  function  of  the  PEDOT  is  chosen  so  that  the  lifetime  vs. 
open-circuit  voltage  predicted  by  Silvaco  is  close  to  that  measured  experimentally.  This  results 
in  a  work  function  difference  of  0.9  Volts.  Unlike  some  work  by  others,  rather  than  injecting 
excitons,  our  simulations  inject  only  free  electron-hole  pairs.  However,  this  seems  unlikely  to  be 
a  major  approximation. 
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5.2  Open-circuit  voltage  simulations  (Pre-radiation) 


To  determine  the  open-circuit  voltage  for  a  given  electron-hole  (e-h)  pair  generation  rate, 
we  set  a  zero-current  boundary  condition  on  one  of  the  contacts  then  calculated  the  resulting 
potential  on  that  contact.  The  zero  current  boundary  condition  is  a  goal  which  the  code  attempts 
to  achieve,  and  is  allowed  to  approach  to  within  a  certain  precision.  The  resulting  current  is 
actually  very  small  rather  than  zero.  We  made  several  calculations  under  these  conditions  and 
incorporating  various  combinations  of  Langevin  and  Shockley-Reed-Hall  (SRH)  recombination, 
including  both  and  neither.  (The  carrier  lifetimes  that  are  calculated  when  both  Langevin  and 
SRH  recombination  are  turned  off  are  due  to  the  field  driving  the  carriers  to  the  contacts  where 
they  recombine,  an  obviously  ideal  situation.)  The  goal  was  to  determine  the  dominant 
recombination  method.  The  most  surprising  result  of  this  simulation  was  that  the  predicted  open- 
circuit  voltage  was  only  slightly  affected  by  the  inclusion  of  Langevin  recombination  as  can  be 
seen  in  Fig. 8.  The  inclusion  of  SRH  recombination  along  with  Langevin  caused  a  displacement 
of  the  curve  and  a  different  slope  similar  to  the  results  seen  by  Mandoc7.  In  comparing  these 
results  to  the  Langevin  formalism  used  by  others  it  is  worth  noting  that  Voc  is  still  linearly 
dependent  on  the  log  of  the  generation  rate.  The  conclusion  of  these  simulations  is  that  for  the 
relatively  thin  samples  upon  which  this  effort  was  based,  neither  SRH  nor  Langevin  is  the 
dominant  recombination  mechanism  when  the  density  of  SRH  recombination  centers  is  low,  that 
the  carrier  recombination  follows  with  the  density  of  SRH  recombination  centers,  and  finally  that 
Langevin  recombination  does  not  play  a  role  in  carrier  recombination. 


5.3  Relaxation  Time  Simulations  (Pre-radiation) 

In  order  to  simulate  the  experimental  extraction  of  e-h  lifetime  we  initially  we  used  a  step 
function  in  which  the  e-h  generation  rate  dropped  from  a  high  value  to  a  low  value.  This  is  in 
contrast  to  the  experiment  in  which  a  constant  light  is  shown  on  the  switch  and  a  pulse  of  light  is 
added  and  removed.  While  this  step  approach  produced  reasonable  results,  we  decided  to 
simulate  the  pulse  nature  of  the  experiment. 
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Figure  8.  Predicted  Voc  as  a  function  of  the  carrier  injection  rate  for  different 
recombination  methods;  (□)  no  SRH  or  Langevin,(>)  Langevin  only,  (T)  Tsrh  =10'6  s,  (A) 

tSrh  =10  7  s,  and  (•)  tSrh  =10'8  s. 

At  a  chosen  e-h  generation  rate,  which  corresponds  to  a  chosen  illumination  level,  we  ran 
the  simulation  until  steady  state  was  achieved.  We  then  applied  a  short  duration  (2  ps)  pulse  and 
extracted  the  lifetime  (x)  from  the  behavior  of  Voc  after  the  pulse  ended.  The  height  of  the  pulse 
was  chosen  so  that  AV0C  caused  by  the  pulse  was  similar  to  the  experimental  value  at  Voc  -0.004 
Volts. 

As  seen  in  Fig.  9,  inclusion  of  Langevin  recombination  again  had  a  small  effect  on  the 
extracted  lifetimes  and  slope  of  ln(x)  vs.  Voc.  Thus,  without  any  SRH  the  slope  was 
approximately  -34.5.  According  to  the  Langevin  analysis  from  Eq.  3,  the  slope  is  -kt/q  =  -38.6 
while  the  experimental  values  depend  on  the  cathode  interface  and  range  from  26.5  to  38.6.  The 
inclusion  of  SRH  recombination  changed  not  only  the  magnitude  of  x  but  also  the  slope  of  the 

o 

t(Voc)  curves,  lowering  them  to  as  little  as  -22.7  when  the  lifetime  was  10'  s.  While  the  traps  in 
these  simulations  were  placed  uniformly  in  the  device,  it  is  possible  to  explain  the  experimental 
change  in  the  slope  of  the  lifetime  curves  as  being  due  to  different  trap  densities  and,  therefore, 
SRH  lifetimes. 
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Figure  9.  The  lifetime  as  a  function  of  Voc  for  different  recombination  methods;  (□)  no  SRH 
or  Langevin,(4)  Langevin  only,  (T)  tSrh  =10'6  s,  (A)  tSrh  =10'7  s,  and  (•)  tSrh  =10'8  s. 


5.4  Charging  Simulations 

In  order  to  test  our  hypothesis  that  the  experimental  irradiation  results,  in  which  Voc 
dropped  as  x  remained  constant,  was  due  to  the  buildup  of  trapped  charge,  we  included  a  small 
region  of  trapped  holes  near  the  simulated  anode.  We  increased  this  trapped  charge  level  to 
simulate  the  effect  of  having  a  greater  total  dose.  The  qualitative  matching  was  good  as  shown 
in  Figures  10  and  11.  In  Figure  11,  the  open  circuit  voltage  follows  the  solid  curve  with  triangles 
from  the  upper  left  to  lower  right  as  the  light  level  is  increased.  In  other  words,  light  level  is  the 
parameter  for  this  curve.  Further,  the  curve  shows  that  as  the  light  induced  open  circuit  voltage 
is  increased,  the  carrier  lifetime  is  reduced.  To  simulate  an  experiment,  a  light  level  is  set,  which 
places  the  start  of  a  simulation  at  a  particular  point  along  this  curve,  and  is  the  initial  point  for  the 
simulation.  A  simulated  radiation  is  applied  by  increasing  the  trapped  hole  density  near  the 
anode,  causing  the  open  circuit  voltage  to  drop  from  the  initial  point  as  shown  in  Figure  11. 
These  simulations  are  shown  in  Figure  1 1  with  open  squares,  open  and  closed  circles.  One  can 
readily  see  that  while  the  open  circuit  voltage  drops,  the  lifetime  as  a  function  of  simulated  dose 
is  flat,  up  to  a  point,  at  which  point  all  three  curves  have  notable  upswings  at  low  biases.  At 
present  it  is  not  clear  if  this  is  an  artifact  of  the  calculation  or  something  real  which  could  be  seen 
experimentally  if  the  total  dose  were  high  enough. 
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From  the  same  simulations  we  were  able  to  link  the  level  of  the  charging  to  the  effect  that 
the  charging  had  on  Voc.  The  three  initial  biases  represented  in  Fig.  1 1  all  have  similar  values  for 
charging  below  3x10  “  cm"  ,  and  begin  to  spread  apart  at  higher  densities.  These  higher  densities 
correspond  to  the  upswing  at  low  Voc. 


Figure  10.  The  simulated  change  in  Voc  due  to  the  presence  of  the  trapped  hole  density;  (□) 
the  Voc  with  zero  charge  is  0.33  V,  (o)  0.40  V,  and  (•)  0.45  V. 


Open  Circuit  Voltage  (V) 


Figure  11.  The  simulated  carrier  relaxation  time  as  a  function  of  Voc:  (A)  the  incident  light 
level  was  adjusted  to  vary  Voc,  (□,  o,  and  •)  increasing  trapped  hole  density  near  the  anode 
with  the  initial  VOC=0.33V,  0.40V,  and  0.45  V  respectively. 


6.0  Conclusions 

This  effort  was  initiated  at  the  request  of  AFOSR  with  the  goal  of  determining  the  effect 
of  radiation  upon  the  performance  of  organic  semiconductor  devices.  The  results  of  this 
experimental  and  modeling  effort  were  then  to  be  transitioned  to  organic  device  investigators  for 
use  in  developing  mitigation  techniques.  The  findings  included  the  discovery  that  while  carrier 


11 

Approved  for  public  release;  distribution  unlimited. 


lifetimes  were  unaffected  by  ionizing  radiation,  degradation  did  occur  and  manifested  itself  in 
the  development  of  an  internal  potential  that  tended  to  reduce  efficiency  of  the  organic 
photovoltaic  devices  tested. 

The  team  recommends  additional  modeling  and  experimental  efforts  to  answer  questions 
that  arose.  Specifically,  was  the  upswing  in  carrier  lifetimes  at  high  total  integrated  dose  levels 
real  or  an  artifact  of  the  simulation.  Additionally,  we  recommend  taking  a  further  step  by 
integrating  “real”  space  conditions  into  our  experiments  and  utilizing  ultra-low  dose  rates  (1 
rad/min  (SiCfi)).  Application  of  realistic  dose  rates  will  enable  a  better  understanding  of  the 
effects  of  radiation  on  the  solar  cells  and  will  be  more  representative  of  what  might  happen  if  the 
cells  are  sent  to  space.  We  also  recommend  investigating  the  device  parameters  that  change  due 
to  radiation.  This  also  includes  a  deeper  look  into  the  changes  in  the  shunt  and  series  resistances, 
and  the  ideality  factor  will  help  understand  the  physics  behind  the  cells. 

There  are  several  more  questions  to  answer,  for  which  the  simulations  will  need  to  be 
improved.  The  experimentally-extracted  lifetime  rises  as  the  open  circuit  voltage  falls  but 
reaches  a  saturation  lifetime  in  the  limit  of  low  initial  Voc.  Currently  our  simulations  have  not 
been  able  to  recreate  this  saturation.  It  is  possible  that  a  different  kind  of  recombination  than 
either  Langevin  or  SRH  is  present.  If  so,  the  lifetime  will  not  depend  on  Voc,  or  by  extension,  the 
carrier  densities  since  the  carrier  densities  depend  on  the  field  and  Voc.  It  is  also  worth  noting 
that  device  simulators  have  some  difficulty  accurately  predicting  results  when  light  intensities 
are  low.  Secondly,  the  sensitivity  of  the  parameters  of  an  organic  solar  cell  with  X-ray 
bombardment  to  the  ratio  of  the  P3HT:PCBM  should  be  explored.  At  present  no  modeling  has 
explained  the  observed  sensitivities.  Further  expansion  of  the  simulations  is  necessary. 
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